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Abstract. In this paper, we improve the previous work on the MHD Alfven wave oscillation model for the neutron 
star (NS) kHz quasi-periodic oscillations (QPOs), and compare the model with the updated twin kHz QPO data. 
For the 17 NS X-ray sources with the simultaneously detected twin kHz QPO frequencies, the stellar mass M 
and radius R constraints are given by means of the derived parameter A in the model, which is associated with 
the averaged mass density of star as (p) = 3M/(47rR"^) ~ 2.4 x 10^^(g/cm"^)(A/0.7)^, and we also compare the 
M-R constraints with the star equations of states. Moreover, we also discuss the theoretical maximum kHz QPO 
frequency and maximum twin peak separation, and some expectations on SAX J1808. 4-3658 are mentioned, such 
as its highest kHz QPO frequency ~ 870 (Hz), which is about 1.4-1.5 times less than those of the other known 
kHz QPO sources. The estimated magnetic fields for both Z sources (about Eddington accretion rate MEdd) and 
Atoll sources (~ VYoMEdd) are approximately ~ 10® G and ~ 10* G respectively. 
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1. Introduction 

The Rossi X-Ray Timing Explorer (RXTE) has observed 
the kilo-Hertz quasi-periodic oscillations (kHz QPOs) in 
the X-ray flux in about 25 accreting neutron stars in low- 
mass X-ray binaries (LMXBs), and the twin peak kHz 
QPOs, the upper and lower frequencies (z^2 and z^i), are 
usually shown in the Fourier power spectrum (e.g., van der 
Klis 2000, 2005, 2006 for the recent reviews). Remarkably, 
the ranges of kHz QPO frequencies are almost homoge- 
neous at about ~200 Hz - 1300 Hz for both the less lumi- 
nous Atoll sources and the bright Z sources (on the defini- 
tion of Atoll and Z, sec Hasinger & van der Klis 1989), and 
they increase with the inferred mass accretion rate, which 
implies that the kHz QPO mechanisms should be related 
to the common properties of both classes of sources (e.g., 
van der Klis 2000, 2005, 2006; Belloni et al. 2005; Zhang et 
al. 2006a). Furthermore, it is also found that the separa- 
tion of twin kHz QPO peaks Av ~ V2 — i'i is not a constant 
and usually decreases with the mass accretion rate (e.g., 
Mendez & van der KHs 1999; van der Klis 2000, 2006). 
However, it is found that the peak separations increase 
with the accretion rate when the kHz QPO frequencies 
are low in Cir X-1 (Boutloukos et al. 2006) and in 4U 
1728-34 (MigHari et al. 2003), which supports the non- 
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linear correlations between the twin kHz QPO peaks (e.g. 
Zhang et al. 2006a) , and it has been implied by the Alfven 
oscillation model by Zhang (2004) and relativistic preces- 
sion model by Stella & Vietri (1999). Moreover, the kHz 
QPO frequencies also follow the rather tight correlations 
between themselves and with the other timing features of 
the X-ray emissions (see, e.g., Zhang et al. 2006a; Psaltis 
et al. 1998; Psaltis et al. 1999ab; Stella et al. 1999; Belloni 
et al. 2002; Titarchuk & Wood 2002), and the ratios be- 
tween twin kHz QPO peaks systematically decrease with 
the frequency with the averaged ratio value of about 1.5 
(Zhang et al. 2006a). 

In order to account for the mechanisms of twin kHz 
QPO phenomena, some viable models have been proposed. 
At the early stage of discovery of kHz QPOs, the sonic- 
point beat-frequency model is proposed (see, e.g.. Miller 
at al. 1998), which predicts a constant Av at the stellar 
spin frequency, however, with the further observations, the 
simple beat model is inadequate for the varied kHz QPO 
separations. Instead, in the developed sonic-point model 
by Lamb & Miller (2001), the authors consider the disk 
flow at the spin-resonant radius to be smooth or clumped 
to interpret why the occurrence of kHz QPO pair separa- 
tion is close to spin frequency as detected in XTE J1807- 
294 (see, e.g. Linares et al. 2005; Zhang et al. 2006b) or 



hall spin irequency as detected m oAA Jl»L)».4-c5bo8 (see, 
e.g. Wijnands et al. 2003) 

Later on, the relativistic precession model is proposed 
by Stella and Vietri (1999), in which the upper kHz QPO 
frequency V2 is identified with the Keplerian frequency of 
an orbit in the disk and the lower kHz QPO frequency Vi 
with the periastron precession of that orbit, and the var- 
ied Ai^ can be consistently explained in this model. More 
recently, Rezania and Samson (2005) propose a model to 
explain the kHz QPOs based on the interaction of accret- 
ing plasma with the neutron star magnetosphere, where 
the matter is accelerated by the gravitational pull of the 
compact object and hits the star magnetosphere with a 
sonic or supersonic speed. Basically the conclusions of this 
model are consistent with the observed data for the ap- 
propriate choices of the parameters. 

Moreover, Titarchuk et al. (1998), Osherovich & 
Titarchuk (1999ab) and Titarchuk & Osherovich (1999) 
have developed the alternative models. They requir(^ the 
lower kHz QPO frequency to be due to the Keplerian fre- 
quency of matter in the disk and the upper kHz QPO fre- 
quency to be the hybrid between the lower kHz QPO fre- 
quency and the stellar spin frequency. Nonetheless, there 
have not yet been any agreements on the origins of the 
QPOs in neutron star and black hole binary X-ray sources, 
nor on what physical parameters determine their frequen- 
cies, which have also been identified with various char- 
acteristic frequencies in the inner accretion flows (e.g. 
Abramowicz et al. 2003ab; Abramowicz 2005; Kluzniak 
et al 2004; Lee et al. 2004; Rebusco 2004; Rebusco & 
Abramowicz 2006; Petri 2005; Horak & Karas 2006). 

This paper is an improved one of the previous work 
on the MHD Alfven wave oscillation mechanism for the 
NS kHz QPOs (Zhang 2004), where we consider the ra- 
dial dependence of the accretion flow velocity, and the 
applications have been done by comparing the model's 
predictions with the updated kHz QPO data. 

The paper is organized as follows: In section 2, the 
overview and scenario of the model are described, and 
the derivations of kHz QPO frequencies are given. The 
further applications of model arc presented in section 3. 
The conclusions and consequences are summarized in the 
final section. As the conventional usage, the Newtonian 
gravitational constant G and the speed of light c=l are 
exploited. 



2. The overview of the improved model 

The previous model for the NS kHz QPOs proposed by 
Zhang (2004) does not mention the excitation mechanism 
of Alfven wave oscillation to modulate the X-ray flux to 
produce the observed kHz QPOs. Therefore, in this pa- 
per we try to answer these questions. To straightforwardly 
grasp our point of views, an imagined geometrical illustra- 
tion of the magnetosphere and accretion disk associated 
with our model is plotted in FIGURE 1. 



z.i. MHU wave and its excitation in the accretion now 

In MHD, the Alfven wave is a transverse wave and prop- 
agates along the magnetic field line, the effects of which 
have been studied in solar physics to interpret the de- 
tected quasi-periodic oscillation of several minutes in coro- 
nal loops (e.g., Roberts 2000; Nakariakov et al. 1999; 
Aschwanden et al. 1999). As known, the coronal loops may 
be set into oscillations with various modes, such as a kink 
mode (Roberts 2000 for a review), thus we assume the 
similar oscillations to occur in the accretion disks with 
the loop length of circumference 27rr at the disk orbit ra- 
dius r. As pointed out, the MHD turbulence by the shear 
flow in the accretion disk (e.g., Ruediger & Pipin 2000) 
will trigger the strong variation of plasma energy density 
and ignite the shear Alfven wave motion along the or- 
bit. In a wave length of circular perimeter, the intensely 
excited shear Alfven wave can dedicate to the observed 
X-ray flux fluctuations, which is then responsible for the 
observed kHz QPO phenomena. 

However, the propagation of the Alfven wave will arise 
its damping because of the dissipation by the X-ray emis- 
sion and the viscous interactions in MHD, which may ac- 
count for why sometimes we only measure a single kHz 
QPO or nothing. From the conventional accretion disk 
prescription, the accretion flow in the equatorial plane of 
disk will drag the polar field lines into the azimuthal di- 
rection, and it is often assumed that the field strength 
of the azimuthal component is comparable to that 
of the unperturbed polar field in the equatorial plane, 
i.e. B^{r) ~ B{r) (e.g., Ghosh & Lamb 1979; Shapiro 
& Teukolsky 1983). 

2.2. On the preferred radius and the coherence 

On the preferred radius r, it is supposed to be a "coher- 
ence" location at where the characteristic Alfven velocity 
VA{r) (frequency) calculated by the plasma mass density 
of quasi-spherical flow coincides with the Keplerian veloc- 
ity V]i{r) (frequency), i.e., VA{r) = V]i{r). We assume that 
at this radius there is a resonance between the Alfven wave 
frequency and Keplerian orbital frequency. 

From the definition, at the Alfven radius Ra, the 
Alfven velocity is the free fall velocity, expressed as (e.g., 
Shapiro & Teukolsky 1983), 

va{Ra) = B{Rk)I^^ttp{Rk) = vs{Ra) , (1) 
where B{r) is the dipole magnetic field and 

p{r) = M/[SrVsir)] , Sr = A^r^ (2) 

is the plasma mass density defined by the spherical ac- 
cretion with the free fall velocity vs_{i) or the Keplerian 
velocity Wk(i'), 

vs{r) = v'2GM/r = y/2vk{r), Uk(r) = ^/GM/r . (3) 



ihereiore, at the preierred radius r = (pKA, where is 
an introduced parameter, the Alfven velocity equals the 
Keplerian velocity, thus, 



VA{r) = B(r)/v/4^p(r) = V20-^/4«k(r) = Vk{r) , (4) 

the 'coherence' radius is defined by the condition of 
V2?i-7/4 = lor(/) = 22/7~1.2. 



2.3. The identifications of twin kHz QPO frequencies 

Therefore, at the preferred radius, the Alfven wave fre- 
quency is 



VA{r) 



VAir) Wk(r) 



27rr 



27rr 



Mr) 



(5) 



which we identify to the upper kHz QPO frequency (see 
also Zhang 2004), or written as 



^2 = I'k = 1850(CA)A:3/2(Hz) 
= 1295(|^)X3/2(Hz), 



(6) 



with the parameters X=R/r and A = {m/RQY^^ where 
i?6 — R/lO^{cm) and m — are the stellar radius R 
and mass M in the units of 10^ (cm) and solar masses, 
respectively. It is noted that the quantity is a mea- 
surement of the averaged mass density of star, expressed 
as 

(p) = 3M/(47rR3) ~ 2.4 x 10i''(g/cm3)(A/0.7)2. 
The reasons for introducing the "non-Keplerian factor" ^ 
in Eq. © come from the facts of the complicated physical 
environments, which do not satisfy the ideal conditions for 
applying the Keplerian frequency, i.e. the point mass or- 
biting around a central gravitational source in a vacuum. 
The modifications to the Keplerian frequency will be taken 
into account if the following situations are considered: (a) 
the influence by rotation or the modification from Kerr 
spacetime to Schwarzschild spacetime; (b) the considera- 
tion of the factual environment around star deviating from 
the ideal situation of a test particle motion in a vacuum; 
(c) the plasma blob moving in a MHD with the strong 
magnetic field. In principle, this "non-Keplerian factor" ^ 
should be not a multiplicative constant. For instance, it 
could depend upon the radius in the disk at which the 
QPO is produced, and hence upon the frequency itself. 
However, we have to stress that it is just the reason of 
simplicity to chose it as a constant. In practice, for the 
mathematical convenience, we still take ^ = 1 to process 
the calculations, and finally the real parameter A can be 
obtained through dividing it by the "non-Keplerian fac- 
tor" ^, the implication of which is discussed in the last 
section. 

Equivalently, the NS radius can be expressed by the 
parameters A and m. 



Re 



1.27mi/3(^/0.7)-2/3 (10 km) 



(7) 



Furthermore, on the formation scenario of the lower kHz 
QPO frequency (i^i), we have the following arguments. 



i he accreted materials accumulated on the polar cap will 
become denser by a factor of ratio of the spherical area 
to the polar cap area. These denser materials will flow 
equator-ward and pass through the perpendicular field 
lines by the instabilities, so this fiow is infiucnccd by the 
Lorentz force. The all denser materials will be expelled 
out through the magnetic tunnel and enter into the tran- 
sitional zone near the 'coherence' radius because the flow 
along the field lines is force free and will not experience the 
Lorentz force. Henceforth, we ascribe the Alfven wave fre- 
quency calculated by the denser mass density of plasma 
formed on the polar cap to be the lower kHz QPO fre- 
quency. For the mathematical convenience, it is presumed 
that both the upper and lower kHz QPOs occur at the 
same radius. Thus, the lower kHz QPO frequency i^i is 
described below, 



VApjr) 
27rr 



(8) 



where the Alfven velocity VAp(r) is defined by the mass 
density expelled from the polar cap, 



VApir) 



with 



Bjr) 




Vkir)^ 



Pp 



Pp = M/[SpVsiR)], 



(9) 



(10) 



where the magnetic polar cap area Sp is obtained as 
(Zhang & Kojima 2006), 



Sp ==47ri?2(i_cos6'), 
with 

sin^e = R/r = X, 



(11) 



(12) 



where 6 is the open angle of the last field line to close 
at radius r. As an approximation, the polar cap area is 
usually written as Sp = if i? <C r (see e.g. Shapiro 

& Teukolsky 1983, P.453). Henceforth, if we ascribe the 
upper and lower kHz QPO frequencies to the two Alfven 
wave frequencies with different mass densities described 
in Eqs.® and then we have 



!1 

1^2 



lp{r) ^ SpVBjR) 
Pp V Srvs{r) 



(13) 



After considering Eq.(IT|) and Eq. (fTT|) with the parameter 
definition X= R/r, Eq.([13l) gives 



(14) 



For convenience, the ratio of the twin kHz QPO frequen- 
cies is obtained to be. 



^ X 



-5/4 



l + y/l-X, 



(15) 



which only depends on the position parameter A= K/r 
where X-ray flux responsible for kHz QPOs emits and is 
independent of the averaged mass density parameter A 
and mass M. Furthermore, the twin kHz QPO separation 
is written as, 



(16) 



which is not a constant with the variation of vi or 1^2 • 

From Eq.® and Eq.lUll), if the twin kHz QPO fre- 
quencies are known simultaneously, then the values of A 
and X can be determined. For the detected sample sources 
listed in TABLE I (see e.g. van der Klis 2000, 2006; Belloni 
et al. 2005; data are provided by T. Belloni, M. Mendez 
and D. Psaltis), such as Sco X-1 (van der Klis et al. 1997; 
Mendez & van der Klis 2000), 4U1608-52 (Mendez et al. 
1998ab) , 4U1735-44 (Ford et al. 1998) and 4U1728-34 
(Mendez & van der Khs 1999). 

The comparisons of the model's conclusions to the well 
detected kHz QPO sample sources are shown in FIGURE 

2, and the agreement between the model and the observed 
QPO data is quite good for the selected ranges of NS 
parameters A=0.6, 0.7 and 0.8. In FIGURE 2b, we find 
that Av increases with V2 if V2 <^ 750 (A/0.7) Hz and Av 
decreases with 1^2 if 1^2 >^ 750 (A/0.7) Hz. In addition, the 
twin kHz QPO ratio (FIGURE 2c) follows the decreasing 
tendency with i>2- 

3. The applications of the model 

To inspect the model's predictions, we demonstrate some 
applications of the model in the following. 

3.1. The constrain conditions of NS mass and radius 

The NS mass constrain condition by the kHz QPOs has 
been given by Miller et al (1998) and Zhang et al (1997) 
through assuming that the accretion disk radius of show- 
ing kHz QPOs is bigger than the innermost stable circular 
orbit (ISCO: three Schwarzschild radii), i.e., the maximum 
observed frequency is presumed to be the Keplerian fre- 
quency at ISCO with the condition that the stellar surface 
is enclosed by ISGO, 



m < 2.2/j/2fc 



(17) 



From Eq.([7]), the new radius constrain condition is ob- 
tained if the mass constrain condition is known, therefore 
we exploit the generally assumed NS mass lower limit 1.0 
Mq from the astrophysical argument and its upper limit 
condition Eq. (jl7p to set the radius constrain conditions, 
respectively. 

Re > 1.27(^/0.7)"2/^ 



and 

i.a.(^)^/^^1.65.-V3(_i)-2/3, 



(18) 



(19) 



and radius relations and their constrain conditions are cal- 
culated for the 17 known sources whos twin kHz QPOs are 
simultaneously detected, which are listed in TABLE I. In 
Figure 3, the mass-radius relations inferred from the val- 
ues of parameter A have been plotted, and it is found 
that, for the conventionally accepted NS mass lower limit 
^ Mq , there exists difficulty in accordance with the many 
modern equations of states (EOSs) if the parameter A is 
too low, for instance A - 0.47 for SAX J1808.4-3658 (see 
TABLE I). In general, the lower the value of A, the more 
difficult it is to reconcile the M-R relations with the re- 
alistic EOSs of stars. The situation can be altered if the 
introduced 'non-Keplerian factor' <^ < 1, which makes the 
parameter A increased. 

3.2. The maximum twin l<Hz QPO frequencies and 
their separations 

The theoretical maximum frequency separation A^'max 
can be calculated from Eq. (fTOI by vanishing the variation 
of Al/ respect to X, c?[Ai/max]/'iA' — 0, where we obtain 
X ~ 0.69, 1/2 = 750(A/0.7)(Hz), vi = 382(A/0.7)(Hz) 
and Aj/jnax = 368 (A/0. 7) (Hz). If 1/2 is lower (higher) than 
750 (A/0.7) (Hz), then Ai/ will increase (decrease) with 
1^2- In the recent observations of 4U 1728-34 (Migliari et 
al 2003), the twin simultaneously detected kHz QPO fre- 
quencies are found at the central frequency vi = 308 
(Hz) and 1^2 = 582 (Hz) ( Av = 274 (Hz)), so this is 
the first detected event for a significant decrease of kHz 
QPO peak separation towards low frequencies, however, 
which is qualitatively consistent with our model's predic- 
tion. Nonetheless, the maximum values of twin kHz QPO 
frequencies coincide and occur at X=l where the preferred 
radius equals the stellar radius in case the ISCO is inside 
the star, i/i^ax = i^2max = 1850A(Hz) = 1295(v4/0.7)(Hz), 
namely 



A> 



1850(Hz) 



(20) 



In Eq. (fT9|) . the radius is constrained by A and P2, or equiv- 
alently by the twin kHz QPOs. Furthermore, the NS mass 



In TABLE I, the theoretical maximum kHz QPO frequen- 
cies of 17 kHz QPO sources are listed, and the averaged 
values of A and the maximum position parameter Xmax 
are given for both the Z and Atoll sources, which are, re- 
spectively, (A)=0.66 and (X,nax)=0.87 (r - 1.15i?) for Z 
sources, and (A)=0.74 and X,„ax=0.9 (r l.li?) for AtoU 
sources. However, we currently cannot figure out what 
physical mechanisms arise these systematical differences 
between the Z and Atoll sources. The inequality X^ax < 1 
tells us the fact that the accretion disk does not reach the 
stellar surface, so it is likely that the ISCO prohibits the 
disk from arriving at the stellar surface. If this scenario is 
plausible, we can conclude that the ISCOs of almost all 
sources in TABLE I locate outside the stars and the esti- 
mated positions of their masses and radii should appear 
above ISCO fine in M-R diagram shown in FIGURE 3. 
Interestingly, the one unusual case is the Atoll source 4U 
0614-^09 (e.g., van Straaten et al 2000; van der Klis 2000, 
2006), and its parameter A ~ 0.76 is inferred from the si- 



muitancously detected twin kHz tj-rUs, which imphcs its 
maximum upper (lower) kHz QPO frequency to be about 
1406 (Hz). However, from the observation, one single kHz 
QPO peak frequency 1330 (Hz) was detected in this source 
(van Straaten et al 2000; van der Khs 2000, 2006), show- 
ing X ^ 0.97, thus this detected single frequency may be 
near the maximum kHz QPO saturation frequency, show- 
ing the X-ray spectrum information near the stellar sur- 
face, the confirmation of which needs the further proposal 
and analysis of this source. 

As for the level-off or frequency saturation of kHz 
QPO, it has been paid much attention since the early 
discovery of kHz QPOs, which is ascribed to the occur- 
rence of ISCO or stellar surface (e.g., Zhang et al 1998; 
Kaaret et al 1999; Miller 2004; Swank 2004). However, 
from our model, if Risco > the saturation frequency 
would occur at ISCO with the maximum upper kHz QPO 
frequency i'2max = 2200/m (Hz) (see also, e.g. Miller et 
al 1998; Miller 2004); if Risco < R, the saturation fre- 
quency would occur at the stellar surface R with i^2max = 
1295(^/0.7) (Hz). Therefore, the model inferred maxi- 
mum positions of showing the kHz QPOs satisfy X^ax < 1 
(see TABLE I) for the simultaneously detected kHz QPO 
sources, which implies that either their ISCOs appear out- 
side stars or the unclear mechanisms diminish the kHz 
QPO X-ray spectra just above the stellar surface, which 
needs the further investigations. 

3.3. The estimations of magnetic field strengths of 
NSs in Z and Atoll sources 

The estimation of the NS magnetic field strength B can 
be given by the definition of its magnetosphere (see, e.g. 
Shapiro & Teukolsky 1983), or described by the accre- 
tion induced magnetic evolution model (see, e.g. Zhang & 
Kojima 2006; Cheng & Zhang 1998), 



R 



(21) 



Bf 



4.3 X 10^ (G) ({M)/MEddf'^m^'^R-'''^ , (22) 



where Rm is the magnetosphere radius defined by the long- 
term accretion rate (M) and MEdd is Eddington limited 
accretion rate. kHz QPOs are assumed to be produced 
around Rm with the variation of the instantaneous accre- 
tion rate, thus i?/RM should be comparable to the aver- 
aged position parameter {X), i.e. {X) ~ i?/RM, which can 
be estimated by Eq.® as 



{^2 



1850A(X)3/2 



(23) 



where {V2) is the averaged value of the detected upper kHz 
QPO of Z (-870 Hz) or Atoll (-980 Hz) sources and both 
sources share the homogeneous kHz QPO distributions 
(see, e.g. Zhang et al. 2006a), therefore. 



B ~ {1850A/{iy2)y/^Bf (24) 



lO^(G) 



So, the magnetic helds are proportionally related to the 
accretion rates, and for both Z sources (Eddington accre- 
tion rate Msdd) and Atoll sources (~ l%MEdd) they are 
about — 10^ G and ~ 10^ G respectively if m ^ 1 and 
i?6 ^ 1, which are consistent with the originally hinted 
values from the X-ray spectra of both sources (Hasinger 
& van der Khs 1989). 

On the correlation between QPO frequency and the 
accretion rate, we can write it in the following by means 
of the definition of magnetosphere (Shapiro & Teukolsky 
1983), 



^ {R/rf/' - {R/Ruf^H^u/rr/' 



.(M) 3/7. 
^ B^ ' ^ 



M . 
{MY 



]3/7 



52 



3/7 



(26) 
(27) 



where ~ M / (Al) is a ratio of the instantaneous accre- 
tion rate to the long-term accretion rate. If the neutron 
star accreted ~ O.OIMq, the NS magnetic field will enter 
into a 'bottom state' where the B-ficld is proportionally 
related to the accretion rate as B oc (M)^/^ \)j |;he accre- 
tion induced magnetic evolution model (Zhang & Kojima 
2006). Thus, we obtain a unified expression of QPO fre- 
quency vs. the accretion rate for both Atoll and Z sources 
to be 1^2 ~ Mx , although both sources share a diversified 
luminosity of even more than two magnitude orders. 



4. Summaries and conclusions 

In the paper, we compare the observations with the im- 
proved model for kHz QPOs, and the main conclusions 
are summarized in the following. 

(1) The theoretical maximum twin kHz QPO frequen- 
cies coincide and occur at about 1295(A/0.7) (Hz) when 
the accreted matters to show these QPOs clash on the 
NS surface, and the maximum twin kHz QPO separation 
is 368 (A/0.7) (Hz). (2) For SAX J1808.4-3658, we find 
that its stellar mass density parameter A=0.47 is about 
1.5 times less than the typical values of other kHz QPO 
sources with A~0.7, and we also obtain its highest kHz 
QPO frequency to be ~870 (Hz), which needs the proof 
of future detections. (3) The averaged mass density of NS 
can be described by the defined parameter A as 
(p) = 3M/(47rR3) ~ 2.4 x 10i'^(g/cm3)(A/0.7)2, and we 
obtain (p) ~ 2.4 x 10"(g/cm3) for most of NS kHz QPO 
sources with A ~ 0.7. (4) With the derived parameter A 
listed in TABLE I for the ideal case of the "non-Keplerian 
factor" f = 1, the mass-radius relation curves are plot- 
ted in FIGURE 3, and we find that the EOSs of strange 
matters (CSl and CS2) seems to be not favorite except 
A ^ 1, resulting in an extremely high kHz QPO frequency 
i^2max 1800 (Hz). Morcover, for the generally assumed 
NS mass lower limit, ^ M© for instance, EOSs of the nor- 
mal neutron matters (CNl and CN2) do not fit for the 
stars in the detected kHz QPO sources unless A ~ 0.88, 
corresponding to J^2max ~ 1600 (Hz). If the EOSs of CPC 
(the star core becomes a Bose-Einstein condensate of pi- 
ous) are the possible choices, the mass and radius ranges 



oi these stars arc troni I.UMq to LyMq and irom 15 km to 
18 km, respectively, for A — 0.45 0.79. In addition, if the 
introduced "non-Keplerian factor" ^ is less than unity, for 
instance ^ ~ 0.7, then the derived values of A in TABLE 
I will be increased by a factor of about 1.4, for instance 
from A ~ 0.7 to A ~ 1.0, which in turn, as shown in 
FIGURE 3, results in a fact that many sources listed in 
TABLE I may be the candidates of strange stars (see, e.g., 
Cheng et al 1998; Lattimcr & Prakash 2004). However, in 
our model it seems to be difficult to take the star in SAX 
J1808.4-3658 (its A=0.47) as a candidate of strange star 
as expected by Li et al. (1999), as seen in FIGURE 3. 
Nevertheless, the EOSs of stars in kHz QPO sources are 
still the open issues before the physical influences of the 
"non-Keplerian factor" arc thoroughly settled. (5) From 
the homogeneous kHz QPO distributions for both Atoll 
and Z sources (see, e.g. Zhang et al. 2006a), we conclude 
the NS magnetosphere scales of both sources to be similar, 
which will arise the NS magnetic field to be proportionally 
related to the long-term averaged accretion rate, i.e. NS 
in Z source possesses a stronger magnetic field than that 
of Atoll source. 

Furthermore, the parameter A = 0.78 is implied by 
the model for KS 1731-260, corresponding to the max- 
imum kHz QPO frequency 1443 (Hz), which is bigger 
than the known detected maximum value U2 ^ 1330 (Hz) 
(4U 0614-1-091, van Straaten et al. 2000). Moreover, the 
proposal for the detection of 1500 (Hz) and 1800 (Hz) 
QPO frequency is suggested by Miller (2004), correspond- 
ing to j4 ~ 0.8 and A ~ 1.0 in our model, respectively, so 
these measurements of QPOs above 1500 (Hz) therefore 
have excellent prospects for stronger constraints on the 
mass and radius relations, as well as on the models. It is 
claimed by Miller (2004) that a QPO frequency as high 
as 1800 (Hz) would be large enough to argue against all 
standard nucleonic or hybrid quark matter EOS, leaving 
only strange stars (see also FIGURE 3). Nevertheless, the 
proposal of detecting the QPO data at low frequency vi 
100 (Hz) is also meaningful, by which the model is tested 
through inspecting Ai' vs. 1^2 relation as shown in the mid- 
dle panel of FIGURE 2. It is also predicted by Stella and 
Vietri (1999) Av will increase with the accretion rate if 
1/2 <~ 700 Hz, and the evidence for this has been recently 
detected in Cir X-1 by Boutloukos et al. (2006). 

As a summary, it is remarked that our model is still a 
simple one, and the further improvements are still needed 
through c;onsidering the details of accretion flow or disk 
structure, where the plausible existence of the transitional 
layer (e.g., Titarchuk et al 1998; Titarchuk & Osherovich 
2000) or the nonlinear disk resonances (e.g., Abramowicz 
et al. 2003ab) may be important. In addition, we have not 
yet considered the rotational effect of Kerr spacetime (see, 
e.g.. Miller 2004; van der Klis 2006), the factual MHD 
flow velocity around star but perhaps not a Keplerian 
velocity of test particle, the accretion flow in a strong 
gravity regime, the plasma instabilities and turbulence in 
strong magnetic field, etc. Therefore, the considerations 



01 above comphcations will improve the present version or 
the model. 
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Fig. 1. An imagined schematic illustration of an accret- 
ing neutron star magnetosphere and disk for the kHz QPO 
productions. is the open angle between the magnetic 
axis and the closed field line of magnetosphere, defined 
by sin^ 9 — R/r, where R and r are the radii of star and 
its magnetosphere, respectively. The quasi-spherical ac- 
cretion from all over star falls in and then cumulates at 
the magnetic polar cap, where the blobs or patches of ac- 
creting plasmas are piled up and condensed into the high 
mass density, and some of materials will diffuse onto the 
whole stellar area through the side-flow on account of the 
plasma instabilities and the others sprout out along the 
closed field lines through the out-flow to enter into the or- 
bital flow. The preferred radius is supposed at where the 
shear Alfen wave frequency with the quasi-spherical accre- 
tion propagating in the orbit equals the orbital Keplerian 
frequency, namely the "coherence" is presumed to occur 
there. 
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TABLE I. Neutron star parameters oi detected kHz QPO sources 



Source* 






m(^) 


R^4) 






X 






m^/^ (km) 


(Mo) 


(km) 


(Hz) 


(Hz) 




Millisecond pulsar 


SAXJ 1808.4-3658" 


0.47 


16.9 


1.0-2.99 


16.8-24.4 


870 


725 


0.90 


Z source 


Sco X-1 


0.66 


13.3 


1.0-2.05 


13.3-17.0 


1221 


1075 


0.93 


GX 340+0 


0.64 


13.5 


1.0-2.62 


13.5-19.0 


1184 


840 


0.80 


GX 349+2 


0.68 


13.0 


1.0-2.23 


13.1-17.2 


1258 


985 


0.86 


GX 5-1 


0.63 


13.8 


1.0-2.49 


13.8-18.8 


1166 


890 


0.84 


GX 17+2 


0.67 


13.2 


1.0-2.02 


13.3-16.8 


1240 


1087 


0.93 


CygX-2 


0.70 


12.8 


1.0-2.19 


12.8-16.7 


1295 


1005 


0.85 


Atoll source 


4U 0614+09 


0.76 


12.2 


1.0-1.65 


12.1-14.4 


1406 


1330 


0.97 


4U 1636-53 


0.74 


12.4 


1.0-1.79 


12.1-15.1 


1369 


1230 


0.94 


4U 1608-52 


0.69 


12.9 


1.0-2.00 


12.8-16.2 


1277 


1099 


0.90 


4U 1702-43 


0.75 


12.3 


1.0-2.03 


12.2-15.6 


1388 


1085 


0.85 


4U 1728-34 


0.77 


12.0 


1.0-1.88 


11.9-14.8 


1425 


1173 


0.88 


KS 1731-260 


0.78 


11.9 


1.0-1.83 


11.9-14.6 


1443 


1205 


0.89 


4U 1735-44 


0.73 


12.5 


1.0-1.91 


12.4-15.5 


1351 


1150 


0.90 


4U 1820-30 


0.73 


12.5 


1.0-2.00 


12.4-15.8 


1351 


1100 


0.88 


4U 1916-053 


0.73 


12.5 


1.0-2.08 


12.4-16.0 


1351 


1058 


0.85 


XTE J2123-058 


0.72 


12.6 


1.0-1.93 


12.5-15.8 


1332 


1140 


0.91 



*: The sources are chosen from those that the twin kHz QPOs are detected simultaneously (van der Klis 2000, 
2006; the original references therein), then 4U 1915-05 is not included because its two incompatible values of 
are reported (van der Klis 2006). a: Wijnands et al. 2003. Calculated by the simultaneously detected twin kHz 
QPO data. Obtained by Eq.lO. Estimated by the generally assumed NS mass lower limit 1.0 and by the 
constrain condition m < 2.2/(^21. (Miller et al. 1998). Obtained by Eq.dlSl) and Eq.jnj. !/2max = 1295(^1/0.7) 
(Hz) is the possible maximum kHz QPO frequency. ' ' : i^2obs is the detected maximum kHz QPO frequency. the 
detection inferred maximum X position, Xmax = (i^2obs/!^2max)^^^. 



